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Abstract. Coherent photoproduction of 7r°-mesons from threshold (Eth ~ 136 MeV) throughout the A- 
. . ' resonance region and of ^-mesons close to the production threshold (E t h ~ 570 MeV for rf) has been 

measured for 7 Li nuclei. The experiment was performed using the tagged-photon beam of the Mainz 
. MAMI accelerator with the Crystal Ball and TAPS detectors combined to give an almost 4tt solid-angle 

electromagnetic calorimeter. The reactions were identified by a combined invariant-mass and missing- 
energy analysis. A comparison of the pion data to plane-wave impulse modelling tests the nuclear mass 
form factor. So far coherent ^-production had been only identified for the lightest nuclear systems ( 2 H 
and 3 He). For 3 He a large enhancement of the cross section above plane- wave approximations had been 
reported, indicating the formation of a quasi-bound state. The present Li-data for 77-production agree with 
a plane-wave approximation. Contrary to 3 He, neither a threshold enhancement of the total cross section 
nor a deviation of the angular distributions from the expected form-factor dependence were observed. 

PACS. 13.60.Le Meson production - 14.20.Gk Baryon resonances with S=0 - 25.20.Lj Photoproduction 
reactions 
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1 Introduction 

Photoproduction of mesons off nuclei involves in general 
many different final states of the meson-nucleus system 

Correspondence to: B. Krusche, Klingelbergstrasse 82, CH- 
4056 Basel, Switzerland, e-mail: Bernd.Krusche@unibas.ch 



and can contribute to a wide range of topics (see [I] for 
an overview). Very interesting for many questions are two 
limiting cases. In 'quasi-free' processes, the reaction in- 
volves one specific nucleon, called 'participant', which is 
kicked out of the nucleus, and the rest of the nucleus can 
be regarded as a 'spectator' system that only compensates 
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the momentum of the bound participant. For light nuclei, 
this process is a powerful tool for the study of reactions 
off quasi- free neutrons [I] ; for heavy nuclei it can be used 
as a testing ground for meson - nucleus interactions and 
hadron in- medium properties [2J. 

In 'coherent' reactions, ideally the meson is produced 
via a superposition of the reaction amplitudes from all 
nucleons and, in the final state, the nucleus remains in 
its ground state. A similar process in which no nucleon 
is removed from the nucleus but the nucleus is excited to 
a higher lying nuclear state, is sometimes called 'incoher- 
ent' production. The advantage of the coherent process is 
the simplicity of the final state; the ground-state proper- 
ties of nuclei are well under control. This reaction is well 
suited for the study of the in-medium properties of mesons 
and nucleon resonances. The undisturbed final state can 
be easily constructed from the plane-wave impulse ap- 
proximation and any deviations may be attributed to nu- 
clear effects like meson-nucleus final-state interactions or 
in-medium modifications of hadron properties. Such pro- 
grams have been pursued in particular for the study of 
medium effects on the production and propagation of the 
Z\-resonance in medium via the coherent jA — > Att° reac- 
tion (see e.g. [3112 )• The same reaction was also exploited 
for the study of nuclear properties such as nuclear mass 
form factors [5] , and in incoherent production, for nuclear 
transition form factors [5]. Nuclear form factors in the 
region of helium and lithium isotopes have gained much 
new interest in connection with the study of halo nuclei 
(see e.g. [HIE])- Coherent pion photoproduction allows the 
direct study of the nuclear mass distribution because pro- 
duction of 7t mesons in the Z\-resonance region couples 
identically to protons and neutrons. 

Until now, coherent photoproduction of heavier mesons 
off nuclei has almost not been investigated since such mea- 
surements are very demanding. For mesons like the 77, large 
momenta are transfered to the nucleus, which suppresses 
the production cross section due to the nuclear form fac- 
tors. Background from breakup reactions, where the par- 
ticipating nucleon is removed from the nucleus, dominates 
the production process. This background must be sup- 
pressed either by detection of the recoil nucleons or by 
conditions on the reaction kinematics, demanding detec- 
tor systems with large solid-angle coverage, large detection 
efficiency, and excellent energy and angular resolution. 

Recently, coherent photoproduction of ?7-mesons from 
light nuclei has attracted interest as a tool for the search 
of so-called 77-mesic nuclei [MTU] . The question is whether 
the strong interaction allows the formation of quasi-bound 
meson-nucleus states, which would be the ideal system 
for the study of meson-nucleus interactions. The interac- 
tion of low-energy pions with nuclei is too weak for quasi- 
bound states but the situation is much different for 77- 
mesons. Production of 77-mesons in the threshold region is 
dominated by the excitation of the s-wave Sn(1535) res- 
onance |11U12) . which couples strongly (branching ratio 
Ps 50% [13]) to Nrj. As a consequence, the interaction of 
?7-mesons with nuclear matter is important also for very 
small momenta of the mesons. Typical absorption cross 



sections are around 30 mb and are over a wide range of 
kinetic energy (r«l MeV - 1 GeV) almost independent 
of T [T4"lll5| . First evidence for an attractive s-wave r]N 
interaction, which might lead to the formation of quasi- 
bound states, was reported from coupled channel analy- 
ses of pion-induced ^-production reactions [T51IT7] in the 
1980s. However, it is still controversially discussed whether 
the interaction is strong enough to form such states. The 
original prediction was for nuclei with mass numbers A 
in the range slightly above 10. However, refined values for 
the jyA^-scattering length extracted from more precise re- 
cent Ty-production data extended the discussion to very 
light nuclei like hydrogen and helium isotopes. (See [TO] 
and refs. therein for a summary of recent results.) 

A much explored experimental approach to identify 
77-mesic states is the study of the threshold behavior of 77- 
production reactions. Quasi-bound states in the vicinity of 
the production threshold should give rise to an enhance- 
ment of the respective cross section over phase-space be- 
havior. Many hadron- induced reactions (see refs. in [10) 1 
have been studied for this purpose. Interesting threshold 
effects have been observed for many of them. Particularly 
strong enhancements were found for the pd — > ?y 3 He [18 
and dp — > ry 3 He reactions [TI?II2(T1[2"T] . implying a large 77 He 
scattering length. If such effects are due to a resonance in 
the 77-nucleus system, they should exist independently on 
the initial state of the reaction. 

Electromagnetic induced reactions, like photoproduc- 
tion of mesons, offer a very clean way to study the 77- 
nucleus final state, but have small production cross sec- 
tions, in particular for the coherent process. Photoproduc- 
tion of 77-mesons in the threshold region has been stud- 
ied for several hydrogen and helium isotopes [9l fTTll22l 
l23 1 l23 P 5 p ) P 7 1 and these results allowed the character- 
ization of the spin and isospin structure of the relevant 
transition amplitudes [28] . The reaction is dominated by 
the excitation of the Sn(1535) resonance via the Eq + - 
multipole, which involves a spin-flip of the participating 
nucleon. This means, that coherent 77-production is prac- 
tically forbidden for nuclei with spin J = ground states. 
Also for nuclei with non-zero ground-state spins, depend- 
ing on the nuclear structure, only a fraction of the nucleons 
(those which can participate in spin-flip transitions) may 
contribute. Furthermore, the electromagnetic excitation of 
the Sn(1535) resonance is mainly isovector (A[^ 2 /A^ 2 s» 

0.1, where A^, 2 is the helicity coupling for the proton and 
Ay 2 is the isoscalar part of the helicity coupling) [28], 
so that contributions from protons and neutrons will can- 
cel to a large extent in coherent 77-production. Together 
with the large momentum transfers involved, these fea- 
tures lead to very small reaction cross sections. 

Only nuclei with ground-state spin J and isospin / 
different from zero are promising candidates for the ob- 
servation of the coherent process. Previous experimental 
results are consistent with this picture. The cross section 
for coherent production off the deuteron (J=l, 1=0) is 
small [24], (typical values for da/dfi are on the order 
of 10 nb/sr). Only upper limits have been extracted for 
the J = I = nucleus 4 He [23]- The most interest- 
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Fig. 1. Total cross section for the 7 3 He — > r; 3 He reaction 
|10j compared to plane-wave impulse approximation. Vertical 
dashed lines indicate coherent and breakup threshold for r\- 
production. Insert: ratio of data and impulse approximation. 

ing case studied so far is the J = / =1/2 nucleus 3 He 
[9,1(3. The coherent process was clearly identified. The 
energy dependence of the total cross section shown in 
Fig- CD 1" is different from the expectation for reaction 
phase-space. A strong threshold enhancement relative to 
the plane- wave impulse approximation (PWIA), similar 
to the results from hadron-induced reactions [TM2TJ1I2T] . 
is observed. The angular distributions close to threshold 
are more isotropic than expected from the shape of the 
nuclear form factor [10]. Both observations together have 
been taken as indication for the formation of a resonant- 
like meson-nucleus state [91HU]. 

So far, this is the only isolated case where coherent r\- 
threshold production off nuclei could be studied. Almost 
nothing is known experimentally about its systcmatics and 
the validity of the simple plane-wave impulse approxima- 
tion used in [9,10]. The present work therefore aimed at 
the measurement of this reaction from a different light nu- 
cleus. Apart from 3 H, the mirror nucleus of 3 He, which, 
however, is difficult to handle as a target, the lightest sta- 
ble isotope with nonzero ground-state spin (J 71 ' = 3/2 _ ) 
and isospin (/ = 1/2) is 7 Li. In the relevant range of mo- 
mentum transfer its squared form factor |30j . which is ex- 
pected to be proportional to the cross section, is roughly 
smaller by an order of magnitude compared to 3 He [25] . 
However, a factor of ~ 3 in counting statistics may be re- 
covered from the target thickness (number of nuclei/cm 2 ), 
making the measurement feasible. 

This paper is organized as follows. The assumptions 
and inputs for the modelling of coherent ir°- and ^-photo- 
production off 7 Li in plane-wave impulse approximation 
are discussed in Sec. [2J The experimental setup is de- 
scribed in Sec. [3] and the data analysis, in particular the 
identification of events from the coherent process, is dis- 
cussed in Sec. |4] The measured cross sections for coherent 
7T° and ^-production are summarized in Sec. \S\ and com- 
pared to the results of the PWIA modelling. 



2 Plane wave impulse approximation 

The PWIA of the coherent meson production follows the 
work of Drechsel et al. [3J , taking into account the specific 
features of the ir° — A and r\ — A final states. The main 
inputs are nuclear form factors and the amplitudes for the 
elementary meson production reactions off the free nu- 
cleon. The elastic charge form factor Fc of 7 Li has been 
measured with electron scattering over a wide range of 
momentum transfer q [301 31 [132"] . Lichtenstadt et al. [32] 
also reported results for the inelastic transition form fac- 
tor Fcx related to the excitation of the 478-keV state in 
7 Li. Since the charge and mass rms radii of 7 Li are similar 
7,33], we can use the charge form factors as basis. How- 
ever, they include the effects from the charge distribution 
of the proton. For the meson production reactions we need 
instead the distribution of point-like nucleons. Therefore, 
the measured charge form factors must be divided by the 
proton dipole form factor F 2 (q 2 ); the ratios are denoted 
by Fc* and Fc x *- Figure [2] summarizes the charge form 
factors and their parametrizations used in the PWIA mod- 
elling. For the elastic form factor, the parametrization of 
Fc* is also shown. The g-dependence of the inelastic form 
factor Fcx for small values of q is approximated by the 
model results cited in [32] . 
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Fig. 2. Elastic charge form factors from Bumiller et al. [31] . 
(black) triangles, Suelzle et al. 30 , (blue) stars, and Lichten- 
stadt et al. (green) squares. Dashed line: parametrization of 
form factor Fc, Solid line: form factor Fc*. (Magenta) dots: 
inelastic form factor for 478-keV excitation 32 . Dotted curve: 
parametrization of inelastic form factor. 

The construction of the transition amplitudes starts 
from the effective total energy W = V s cS of the inci- 
dent photon (four- momentum P 7 , laboratory energy P 7 ) 
and an off-shell nucleon (four- momentum P/v) with three- 
momentum pn from its motion inside the nucleus 

s cff = (P 7 +Pv) 2 . (1) 

The nucleon momentum pjy is obtained in the factoriza- 
tion approximation [3] from the momentum transfer q to 
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the nucleus by 



Pat 



A-l 



(2) 



where A is the nuclear mass number and all momenta are 
in the laboratory frame (note that the expressions in [3] re- 
fer to the center-of-momentum frame) . The amplitudes of 
the elementary reactions are then evaluated at W(E^, q). 

The amplitude for meson photoproduction off nuclei is 
in general given by 



when all multipoles except the leading M\ + are neglected. 
The incoherent excitation of the nucleus is included, but 
the contribution turns out to be negligible (see Sec. I5.1[) . 
The amplitudes M\ + are again taken from the MAID- 
model [55] . For the full PWIA cross section the incoherent 
sum 
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T = L + icrK, da = \L\ 
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(3) 



with the spin-independent part L and the spin-dependent 
part K . It is efficiently evaluated in the CGLN param- 
eterization [34], involving the four invariant amplitudes 
Fi,...,Fi. 

The simplest case is coherent 7r°-production from spin 
J = nuclei in the A- resonance region [3] . The elementary 
production amplitudes are identical for protons and neu- 
trons. The dominant contribution to — > ir°A for spin 
J = nuclei thus involves the spin/isospin- independent 
part of the production amplitude. In the CGLN repre- 
sentation a spin-independent piece arises only from the 
term with the F% amplitude. Due to the pseudoscalar na- 
ture of the pion and the overall symmetry of the prob- 
lem, this term has a sin((9*) factor (<9*: pion polar angle 
in photon- nucleus cm-system) in the amplitude |5J. Since 
the dominant excitation of the Z\-resonance is not isospin 
dependent, all amplitudes from protons and neutrons add 
coherently, which is reflected in a factor A in the ampli- 
tude. The full evaluation of the L-piece gives: 

da _lq* 



is used. 

The situation is different for 77-production. Since the el- 
ementary reaction is completely dominated by an isovec- 
tor, spin-flip amplitude, there is no piece corresponding 
to Eq. 2] in pion production. Like in the 3 He case |10) . 
the main contribution to coherent production comes from 
the Sn excitation of the unpaired nucleon via a spin-flip 
transition. The main difference is that 3 He has an un- 
paired neutron while 7 Li has an unpaired proton. We use 
therefore a similar PWIA as for 3 He in [TP] , including the 
incoherent excitation via the Fc x *-teiui 



da 



dn 



'(A) AN)' 

k\ ^ qi) \ 



dn 



{FlM 2 ) + F* Cx M 2 )) (7) 



with a parameterization of the measured "fp — » pr\ cross 
section from [36] for the elementary cross section da e \ em . 
The change of phase space between the different cm. sys- 
tems is derived from the photon and 77 three-momenta 



in the photon-nucleon (fc 7 
q[f^) cm. systems. 



W (N) 



q\, ), and photon- nucleus 



7777 = o^\F2 (W)\ 2 A 2 sin\e:)F^) , (4) 3 Experimental setup 



dn 2 k, 

where the ratio of pion and photon momenta q*, k* gives 
the phase-space factor for the photon-nucleus system. Nu- 
merical values for the CGLN amplitude F2 were taken 
from the MAID analysis of pion photoproduction [55]. 

The 7 Li case is complicated by the unpaired proton 
in the lps/2 orbit, which gives rise to additional contri- 
butions involving also spin-flip amplitudes that may con- 
tribute to all four CGLN amplitudes. Apart from elastic 
reactions, the IP3/2 proton may be excited to the lpi/2 
orbit, populating the low lying l/2~ state of 7 Li with an 
excitation energy of 478 keV. Incoherent pion production 
to this final state cannot be separated experimentally from 
the coherent process and is thus included in the measured 
cross sections. The spin-dependent contribution must be 
small compared to the spin-independent contribution be- 
cause it is lacking the A 2 factor, but it is important for 
extreme forward or backward angles (because it has a 
cos 2 (<9*) dependence instead of the sin 2 (6>*) for the spin- 
independent part). These contributions are approximated 
from the leading M\+ multipole. Evaluation of the mul- 
tipole expansion of the CGLN amplitudes for the spin- 
dependent part of the cross section leads to 



da 



im - |j|M 1+ (^)| 2 COS 2 (0:) (F^(q 2 ) +F^(q 2 )) 



(5) 



The experimental setup was identical to the one used in 
[37U38j , apart from the target (liquid hydrogen targets for 
[37U38j . threshold settings, and trigger conditions. 

The measurement used the tagged photon beam [411 
[42] from a primary 883 MeV electron beam of the Mainz 
MAMI accelerator [551110] . The photons irradiated a 7 Li 
target (enrichment 99%) of 5.4 cm length and a density of 
0.534 g/cm 3 , corresponding to a surface density of 0.264 
nuclei/barn. The reaction products were detected with an 
electromagnetic calorimeter composed of the Crystal Ball 
(CB) 05J and TAPS detectors pjl35]. The 672 Nal crys- 
tals of the CB covered the full azimuthal angle for polar 
angles between 20° and 160° around the target, which was 
mounted in the center of the CB. TAPS covered polar an- 
gles between 1° and 20° as a hexagonal wall of 510 BaF2 
crystals, mounted 1.75 m downstream from the target. In- 
dividual plastic detectors in front of each crystal were used 
for charged particle identification. A schematic view of the 
setup, which covered ps 98% of 47r, is shown in Fig. [5J It 
was complemented by a cylindrical Particle Identification 
Detector (PID) JUj] , mounted around the target inside the 
CB, which covered the same solid angle as the CB. 

The experiment trigger was based on a subdivision of 
the CB and TAPS into logical sectors. For TAPS these 
were eight sectors of 64 modules in a pizza-like geometry, 
and for the CB 45 rectangles. The trigger required signals 
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Fig. 3. Setup of the electromagnetic calorimeter combining 
the Crystal Ball and TAPS detectors. Detectors for charged 
particle identification were mounted in the Crystal Ball (PID 
and MWPC) and in front of the TAPS forward wall (TAPS 
Veto-detector) . 



in at least two logical sectors of the calorimeter above a 
threshold of 20 MeV and an analog energy sum of the CB 
modules above 50 MeV. Once a valid trigger had been gen- 
erated, thresholds for the readout of individual modules 
were 5 MeV in TAPS and 2 MeV in the CB. 



4 Data analysis 

The different analysis steps for the identification of pho- 
tons, charged pions, and recoil nucleons are discussed in 
more detail in [571138) . The analysis of coherent neutral 
meson production off nuclei is special in so far as no charged 
particles (no charged pions, no recoil protons) may occur 
in the final state. Detection of charged particles was only 
used to veto events, which simplifies the analysis (there 
was no need to separate charged pions from protons or to 
extract energy information for the charged particles). Ac- 
cepted were only events with exactly two photons (from 
the 7T° — > 77 or rj — > 77 decays) or with exactly six pho- 
tons (77 — > 37T° — > 67). These are particularly clean data 
samples. 

The invariant-mass spectrum of photon pairs for inci- 
dent photon energies below 300 MeV is shown in Fig. |4l 
It is practically background free. No other reactions with 
significant cross section produce two or more photons in 
this energy range. 
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Fig. 4. Invariant mass spectrum for two-photon events for inci- 
dent photon energies below 300 MeV. Statistical uncertainties 
smaller than symbol sizes. The solid (red) curve is a Monte 
Carlo simulation of the detector response. 
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Fig. 5. Missing energy analysis for single 7r° production for 
different incident photon energies. Black dots: measurement 
(statistical uncertainties smaller than symbol size), solid (red) 
curves: MC-simulation for coherent events, dashed (green) 
curves: MC for breakup events, dotted (blue): sum of both. 
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Fig. 6. Invariant-mass spectra for two-photon (left-hand side) 
and six-photon (right-hand side) events in the energy region 
of the ^-production threshold. Solid (red) lines: signal shapes, 
dashed (green) lines: fitted background, dotted (blue) curves: 
sum of both. Ranges of incident photon energies are given at 
left-hand side. 



Double 7T° production sets in with a very low cross 
section around 300 MeV and loss of two of the four decay 
photons is unlikely. The only possible background source 
is production of single 7r° off quasi-free neutrons with loss 
of one decay photon and misidcntification of the neutron 
as photon. However, the corresponding recoil neutrons are 
mostly emitted to forward angles and can be identified in 
TAPS with time-of-flight versus energy and pulse-shape 
analyses. The important step is then the separation of the 
coherent reaction from breakup reactions with emission 
of recoil nucleons. The suppression of such events by the 
required non-detection of recoil nucleons is limited since 
the detection efficiency for recoil neutrons is only on the 
order of 30% (larger than 90% for recoil protons). We use 
therefore in addition the overdetermination of the reaction 
kinematics of the two-body final state. The laboratory ki- 
netic energy of the meson £^ b is directly measured with 
the calorimeter, and its kinetic cm-energy E^ follows from 
the incident photon energy £L. 
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Fig. 7. Missing-energy spectra for events in the r\ invariant 
mass peaks for different ranges of incident photon energy. No- 
tation for curves is as in Fig. [4] Vertical dotted lines: expected 
positions of coherent peaks. Left-hand side: two-photon events, 
right-hand side: six-photon events. 



The mesons are boosted into the cm-system and the 
difference AE n of the two kinetic energies in the cm- 
system is constructed as 



AE m - E* n (E£ ) - E* m (E. 



(8) 
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The result of this analysis for pion production is shown 
in Fig. [5] The peaks at zero missing energy correspond 
to coherent production and dominate the process at low 
incident photon energies. At higher incident photon en- 
ergies breakup background appears at negative missing 
energies. The shape of the signals was generated with a 
full Monte Carlo simulation of the experiment using the 
GEANT3 package 47 . The event generator for the coher- 
ent process was based on trivial two-body kinematics; for 
the breakup reaction, the momentum distribution of the 
bound nucleons was taken into account. Final state inter- 
actions were not taken into account, which explains the 
deviations between data and Monte Carlo in the tails of 
the distributions for higher incident photon energies. The 
separation of coherent and breakup processes, which must 
be done in dependence on the pion angles, is straightfor- 
ward for the energy range up to E 1 — 300 MeV as shown 
in Fig. [SJ At higher incident photon energies, the contri- 
bution from breakup reactions becomes dominant in the 
angle-integrated missing-energy spectra. Across the angu- 
lar distribution missing energy spectra vary. The fraction 
of coherent events compared to breakup is larger for for- 
ward angles, but due to kinematics the separation between 
coherent and breakup events in missing energy is better 
at backward angles. The 'coherent reaction' includes in- 
coherent excitation of the 478-keV level, which cannot be 
resolved by the missing energy analysis. 

The analysis for coherent ^-production follows the same 
scheme. Invariant-mass and missing-energy spectra are sum- 
marized in Figs. [5] and [7] The main difference to tt° pro- 
duction is that, near threshold, the ratio of coherent to 
breakup cross sections is much less favorable. This comes 
from two effects discussed in Sec. [2] The involved momen- 
tum transfers are much larger, suppressing the coherent 
cross section via the form factor. Since furthermore (apart 
from small components in the nuclear wave functions) only 
the IP3/2 proton contributes, the A 2 factor is missing in 
the coherent cross section. 

The invariant-mass peaks from the two-photon decays 
show some background (double it production with two 
undetected photons, single tt° production off quasi- free 
neutrons with one undetected photon and a misidenti- 
fied neutron), which must be subtracted. The invariant- 
mass signals of the six-photon decays are much cleaner. 
In this case, the invariant masses of the three 7r°-mesons 
are also used to identify the reaction as discussed in [TU] . 
The contribution of breakup background to the missing- 
energy spectra is substantial. A clean coherent signal ap- 
pears only in the immediate threshold region. At higher 
energies the signal can be extracted only by fitting the 
simulated line-shapes to the data, which for incident pho- 
ton energies above 650 MeV becomes unfeasible. However, 
the simultaneous extraction of the cross section from the 
two different 77-decay channels gives some estimate for the 
typical level of uncertainty. 

Absolute cross sections were extracted from the mea- 
sured yields with the target surface density, the incident 
photon flux, and the simulated detection efficiencies. The 
latter were generated with GEANT3 {47, simulations. Typ- 



ical values (depending on incident photon energy and po- 
lar angle of the meson) are 20% - 50% for coherent ir° 
production and 35% - 40% for coherent 77-production to 
the six-photon final state and 60% - 70% for the two- 
photon final state. The uncertainty for the detection ef- 
ficiency simulations is smaller than in 38 for two rea- 
sons. Only photons had to be detected, for which the re- 
sponse of the detector system is best understood. There 
is no additional uncertainty from the properties of the 
event generator because in both cases only trivial two- 
body kinematics is involved in the final state. We estimate 
the systematic uncertainty of the detection efficiency be- 
low the 5% level. The incident photon flux was determined 
from the counting of the number of deflected electrons in 
the focal plane by live-time gated scalers. The fraction 
of correlated photons that pass the collimator and reach 
the target (tagging efficiency, ~ 50% for this experiment) 
was determined with special experimental runs. A total 
absorbing lead-glass counter was moved into the photon 
beam at reduced intensity of the primary electron beam. 
The intensity was reduced at the electron source, so that 
no accelerator parameters differed from normal running. 
In addition to these periodical absolute measurements the 
intensity was monitored in arbitrary units during normal 
data taking with an inonization chamber at the end of 
the photon-beam line. The systematic uncertainty for the 
flux measurement is estimated below the 5% level. The 
systematic uncertainty of the surface density of the solid 
7 Li target is estimated as 3% (due to a somewhat irregular 
shape of the target). 

The largest uncertainty is related to the separation of 
coherent signal and breakup background. For coherent 7r° 
production we estimate a systematic uncertainty due to 
this effect of 2% - 5% for incident photon energies from 
threshold to 200 MeV, 5% - 8% between 200 MeV and 
300 MeV, and 8% - 20% between 300 MeV and 500 MeV. 
For 77-production most of this uncertainty is reflected in 
the statistical uncertainties of the yields, which include 
the uncertainty related to the fitting of the missing energy 
spectra. 



5 Results 

The results for the two reaction channels are of different 
quality and intended for different purposes. The 7r°-data 
have excellent statistical quality. In most figures their sta- 
tistical error bars are smaller than the symbol sizes. Al- 
though we compare them here only to PWIA approxima- 
tions to discuss their most important features, they may 
serve as precision tests for more advanced models, taking 
into account the correct nuclear structure of Li and the 
nuclear effects beyond PWIA. 

The pioneering results for 77-production, at a cross sec- 
tion level of 10 - 20 nb, have limited statistical precision, 
but still allow a comparison of the threshold behavior to 
the 3 He case. 
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Fig. 8. Angular distributions for coherent 7r°-production for different ranges of incident photon energy. Curves: results of PW1A 
model normalized in absolute scale to experiment. 
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5.1 Coherent 7r°-photoproduction 

Angular distributions and the total cross section for the 
7+ 7 Li— > 7 Li + 7T° reaction are summarized in Figs. |8j and 
O We discuss first the total cross section. The energy de- 
pendence and absolute magnitude reflect the properties of 
the elementary production cross section off the nucleon, 
trivial factors like A 2 and sin 2 ((9*), the nuclear form fac- 
tor, FSI effects, and possible in-medium modifications of 
the involved nucleon resonances (here the Z\(1232)). We 
compare the data to similar results for the deuteron [15] 
and 12 C 0]. The systematic evolution of the A- resonance 
peak in dependence of the nuclear mass number from 'al- 
most free' production for the deuteron to 'almost nuclear 
density' for carbon is clearly visible. 




200 300 400 500 



E y [MeV] 

Fig. 9. Total cross section for coherent 7r°-production. The 
shaded (green) band indicates the size of systematic uncer- 
tainty of the data. Data for the deuteron (scaled down by fac- 
tor of two) [H] and 12 C [I] for comparison. Solid curve: PWIA 
results, Eq.[6] dashed curve: predicted contribution of spin-flip 
amplitude (first term of Eq. [5]) scaled up by a factor of 10, dot- 
ted curve: incoherent contribution from excitation of 478-keV 
level (second term of Eq. O scaled up by factor of 100. Insert: 
ratio of measured cross section and PWIA, dashed lines: range 
of systematic uncertainty of 7 Li data. 

One should keep in mind, as discussed in detail in [4], 
that the effective position of the Z\-resonance peak is de- 
termined by different effects: the interplay between the 



nuclear form factor and the sin 2 ((9) term in the PWIA 
approximation (see Eq. |4] not valid for the J = 1 deuteron 
which is lacking the sin 2 term), the FSI effects in distorted- 
wave impulse approximation (DWIA), and the density de- 
pendent in-medium modification of the position and width 
of the resonance. Actually, the model of Drechsel et al. [3J, 
which reproduced quite well the data for nuclei from car- 
bon to lead [I], predicts an upward shift of the Zi(1232) 
in-medium resonance position; although the peak in the 
cross section appears to be downward shifted due to the 
other effects. The lithium case is interesting because it is 
transitional between the Z\-in- vacuum and Z\-in-normally- 
dense-matter cases. Previous results [4] have shown that 
the measured cross sections from carbon to lead can be 
reproduced with Z\-self energies extracted from 4 He data. 
However, 4 He is itself a very dense nucleus and the effec- 
tive density of Li is significantly lower than for any of 
the nuclei studied so far. The extraction of Z\-self energies 
from the lithium data will require detailed model calcula- 
tions, taking into account the FSI effects, which are not 
yet available but in progress. 

Here, we compare the measured cross sections to the 
PWIA modeling discussed in Sec. [2] It is obvious from 
the figure that the elastic spin-flip-term (term with Fq^ in 
Eq. [5]) and the incoherent excitation of the 478-keV state 
of the 7 Li nucleus (term with F^ xt in Eq. are negligible 
effects for the total cross section, both much smaller than 
the systematic uncertainty of the data. (Note, however, 
the importance of the spin-flip-term for the angular dis- 
tributions discussed below.) In the low-energy range, up 
to incident photon energies of ps 225 MeV, the measured 
cross sections agree surprisingly well with the PWIA re- 
sults (mostly within systematic uncertainties of the data) . 
This demonstrates that the trivial effects of the coherent 
process are well understood in PWIA and that in this 
regime effects from FSI and in-medium modifications of 
the Z\-resonance must be either both small or cancelling. 
In the maximum of the Zi-resonance, PWIA largely over- 
estimates the data. This is consistent with the expected 
onset of strong FSI and the in-medium damping of the 
Z\-resonance. At even higher incident photon energies, be- 
yond the energy range where the elementary cross section 
is dominated by the Z\(1232) excitation, the model is miss- 
ing contributions from other photoproduction multipoles 
(e.g. from the excitation of the Pn(1440) and Di 3 (1520) 
resonances and background terms), so that no agreement 
can be expected. 

The shape of the angular distributions in Fig.|5Jis quite 
well reproduced at low incident photon energies and even 
reasonably well at higher energies. This is so, because the 
shape is dominated by the sin 2 ((9) term and the nuclear 
form factor. However, a closer inspection of the angular 
distributions also shows some systematic deviations be- 
tween experiment and PWIA for the energy range where 
FSI effects seem to be small. For a more detailed analysis 
Figure [TU] shows a reduced version of the differential cross 
sections as a function of the squared momentum transfer 
q 2 . The cross sections have been divided by the PWIA 
estimate from Eq. [?J] but without the form-factor terms 
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Fig. 10. Form factor of 7 Li extracted from the ratio of mea- 
sured angular distributions and PWIA results for different 
ranges of incident photon energy (see text). The absolute scale 
corresponds to the 152 MeV data, the other data are scaled 
down by successive factors of two. The solid lines correspond 
to fits with Eq. 1101 The insert shows the rms mass radii (red 
dots) extracted from the fits with Eq. 1121 The solid line rep- 
resents the average (dotted lines statistical uncertainty), the 
dashed line the rms charge radius (for point-like protons). 

in Eqs. 2] and [5] The square roots of these ratios, shown 
in the figure, correspond to the nuclear mass form factor 
when the PWIA is valid (and the incoherent excitation 
can be neglected). Shown are only the results for pion cm- 
angles with cos((9*) > —0.5, where the PWIA approach 
seems to be reasonable. The first important observation 
is that the g 2 -dependence of these distributions is almost 
independent of incident photon energy. This is what one 
would expect for a q 2 dependence related to the nuclear 
form factor. 

It was then tested whether the data can be fitted with 
a model of the form factor. The form corresponding to a 
simple harmonic oscillator shell model 



F H o(q 2 ) = d(l - cq 2 )eM-aq 2 ) 



(9) 



did not give satisfying results for the whole range of mo- 
mentum transfers (see below). Much better results were 
obtained with the double- well form for s- and p-orbits used 



Fmo — a o 



2 , 
3 eXp( " 



■9^/4) 



+ g(l - g'a2/6)exp(-^/4) 
b 2 = a 2 {l-l/A) 1 i = l,2 



(10) 



where ciq accounts for the overall normalization and a\, 
ai are the well-strength parameters of the s- and p- wells. 
Fits with this model form factor are shown in Fig. [TTJ 
as solid lines. They excellently describe the data over a 
large range of incident photon energies and momentum 
transfers with almost identical parameters. The average 
values of the well-strength parameters are 



01 = (1.599 ±0.001) fm 
a 2 = (2.47 ±0.06) fm . 



(11) 



Suelzle et al. (3D] quote for the charge distribution pa- 
rameters ai=(1.55±0.015) fm and a 2 =(2.02±0.06) fm so 
that the s-well strength is very similar for charge and mass 
distribution (3% difference), while the p-well strength is 
« 20% larger for the mass distribution. The rms radius is 
related to these parameters by: 



A - I 
A 



1 



1 



(12) 



The values for r rms obtained from the fits are shown in the 
insert of Fig. 1101 They show no systematic variation with 
incident photon energy and their average of ~ 2.62 fm 
is significantly larger than the r rms radius of the charge 
distribution (« 2.27 fm in (3D]; note that the value of 
2.43 fm quoted in this reference includes the charge radius 
of the proton). 

For a more detailed analysis the average of the distri- 
butions from Fig. 1101 for incident photon energies up to 
280 MeV (after renormalization of their absolute scales) 
is compared in Fig. [TT] to the charge form-factor values 
from 30,31,32 . It is evident that the g 2 -dependence of 
the electron scattering data is different from the present 
results. Both data sets have been fitted for the range of 
q 2 < 3 mi" 2 with Eqs. [9] and US The fits with the simple 
harmonic oscillator model (Eq. [5]) were of much inferior 
quality (reduced \ 2 sa 880 for present data compared to 
8 for the double- well form Eq. 1101) and were not further 
considered. The fits with the double- well form from Eq.fTUl 
are shown in Fig.[TTJas solid blue (present data) and solid 
black (electron scattering data) lines. They correspond to 



the following rms-radii (riml '■ electron data 
data): 



rms 



(2.30 ±0.02) fm, 
(2.66 ±0.01) fm. 



present 

(13) 
(14) 



The insert of Fig. [TTJ shows the ratio of the present 
data and this fit (filled, red points). For g 2 -values up to 
3 fm -2 the fit reproduces the shape to within ±0.5%. Due 
to this small systematic differences between fit curve and 
data, the result for the radius is almost independent on 
the fitted range. If, for example, we fit only the data for 
incident photon energies below 225 MeV, where agreement 
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Table 1. Fit results for the mass rms-radius. Results are given in column (1) for fits with the full PWIA model fEqs. I4l5l6)l . 
in column (2) for a truncated model without the spin-flip contribution (Eq. [5} , and in column (3) for a model with the spin-flip 
contribution arbitrarily doubled. Column (4) shows for comparison results of fits to charge form factor from electron scattering 
(divided by proton charge form factor). First row 1 ' average of the fit results from Fig. [TO] with the double- well model (Eg. 1 10p 
over an range of incident photon energies from 150 MeV - 360 MeV (x 2 values are averages for all fits). Second row ' fit with 
the double- well model to the averaged form factor for q 2 < 3 fm~ 2 fFig. lll[) . Third row 3 ' fit with series (Eq. [16] with N = 5 (for 
the model with neglected spin-flip term with N = 7, since N — 5 did not converge). Fourth row 4 -* fit with series with N = 2 
for q 2 < 0.5 fm -2 (only few data for charge form factor). 





full model 


no spin-flip 


spin-flip doubled 


charge form factor 


method 


rims [fm] 


x 2 


Trms [fm] 


x 2 


rims [fm] 


x 2 


riml [fm] 


x 2 


^ double well 


2.618±0.004 


3.4 


2.710±0.004 


7.2 


2.587±0.004 


3.2 






2 ' double well 


2.659±0.007 


8.6 


2.898± 0.003 


24 


2.612±0.002 


8.4 


2.30±0.02 


3.7 


3) series N=5 (7) 


2.635±0.002 


8.4 


2.981± 0.002 


16 


2.575±0.002 


8.3 


2.17±0.04 


1.9 


4) series N=2 {q 2 < 0.5 fm -2 ) 


2.56±0.12 


5.9 


3.12±0.08 


14 


2.398±0.15 


7.2 


2.2±1.2 


1.2 



present exp. 




10 ' 1 

q 2 [fm" 2 ] 

Fig. 11. Average of the form factors from Fig. 1101 compared 
to the charge form factors from electron scattering (divided by 
proton charge form factor) . Solid (blue) line: fit to present data 
(q 2 < 3 fm" 2 ) with Eq.[TDl Results for fit with Eq.QHand N = 
5 (q 2 < 3 fm -2 ), and fit with N = 3 for q 2 <0.5 fm~ 2 (dashed 
and dotted lines) are not distinguishable from solid line; (black) 
solid line: fit of electron scattering data with Eq. 1101 



between data and PWIA is best, the radius changes only 
from 2.659 fm to 2.653 fm. Also shown in the insert (black, 
open points) is the result from an analysis that neglected 



the spin-flip term (Eq. [5]) in the elementary production 
cross section. The influence of this term is substantial at 
small q 2 values; the reduced x 2 of the fit rises from 8.6 to 
24 if it is omitted. 

The rms radius can be also extracted from the present 
data without the use of a specific model for the form factor 
from its slope for q 2 — > 0, using the expansion 



F(q 2 



1 



2 

T 2 

— r 

g rms 



0(q 4 



The data were fitted with the ansatz 

N 



F(q 2 )^J2 c ^ 2n 



71=0 



from which the rms radius follows as 



= y/-6d/ 



CO 



(15) 



(16) 



(17) 



where for correctly normalized form factors Co would be 
unity (here it differs by a few per cent from unity). 

Different fits have been exploited. Two extreme cases 
are fits for the q 2 range up to 3 fm~ 2 with N = 5 and 
with N = 2 only for small momentum transfers (q 2 < 

0.5 fm -2 ). The results for rfml extracted from Ea.[T7lare 
in agreement and close to the above value from the double- 
well harmonic-oscillator model: 



(2.635 ± 0.002) fm 



for the N = 5 fit over the full range and 



(2.56 ±0.12) fm 



(18) 



(19) 



for the slope from the low-momentum transfer N = 2 fit 
(quoted uncertainties are statistical). The fit curves are 
so similar to the double- well result that they are indistin- 
guishable from it in Fig. [TT] 

The results for all model fits are summarized in Ta- 
bic [TJ The form factors derived in PWIA from the coherent 
pion data correspond to an rms-mass radius of (2.60 - 
2.65) fm~ 2 (column (1) of the table), which is significantly 
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larger than the result for the charge radius (column (4) of 
the table). The reduced \ 2 of the fits is larger than unity, 
which is due to the systematic structure of the form factor 
at the sub-percent level (see insert of Fig. [TT1) . which is sig- 
nificant within statistical uncertainties, but much smaller 
than (energy dependent) systematic uncertainties not in- 
cluded in the fitting process. 

One source of systematic uncertainty is the contribu- 
tion of the spin-flip-term (Eq. [5J in the PWIA approx- 
imation, which includes only the dominant M\ + ampli- 
tude and ignores all other production multipoles. Its in- 
fluence has been tested by model fits: excluding it com- 
pletely (column (2) of Table [J) and arbitrary doubling its 
strength (column (3)). Excluding the spin-flip contribu- 
tion increases significantly the \ 2 of the fits and increases 
the value of the radius. Enhancing the spin-flip term by a 
factor of two is quite a large (probably unrealistic) varia- 
tion, since in the Zi-resonance range it is strongly domi- 
nated by the well-known M\ + multipole. The y^-vahxes of 
these fits are similar to the standard version and the radius 
becomes smaller, but is still larger than the charge radius. 
The comparison gives some indication of the possible size 
of systematic uncertainty due to this term. 

The contribution of the inelastic 478-keV excitation 
was ignored for the form-factor extraction. To test its im- 
portance we subtracted the PWIA estimate for this pro- 
cess from the measured angular distributions and repeated 
the analysis. This removes strength at large q 2 , which 
makes the form factor steeper and thus tends to increase 
the radius. However, the effect is smaller than statistical 
uncertainties and can be safely neglected. 

So far no model results are available for FSI effects in 
7 Li. Results for other light nuclei ( 4 He, 12 C) [li] have 
shown that they are important for the energy dependence 
of the total cross section. Nevertheless, the good agree- 
ment of the measured total cross section with the PWIA 
modelling at incident photon energies below 225 MeV indi- 
cates that they must be small for 7 Li in this energy range. 
The main FSI effect depends on the pion kinetic cm ener- 
gies (and thus on the incident photon energy) but it could 
also modify to some extent the shape of the angular dis- 
tributions, which are the basis for the form-factor extrac- 
tion. However, the form-factor fits (see insert of Fig. HTJ1) 

give consistent results for the mass radius r^ml over a wide 
range of incident photon energy, over which the energy-de- 
pendent FSI effects change drastically, from a few per cent 
between 180 MeV and 220 MeV to almost 40% around 280 
MeV (see insert of Fig. 

In order to explain the observed difference between 
the extracted form factor and the charge form-factor data, 
FSI effects with a very peculiar behavior would be needed. 
This is demonstrated in Fig. 1121 where the low-energy an- 
gular distributions are compared to a modified PWIA. 
The only difference to the PWIA curves in Fig. [5] is that 
instead of the form factor from electron scattering the 
double-well parameterization of the present form factor 
data from Fig. [TTJ corresponding to r™l=2.66 fm was 
used. The absolute scales of the PWIA results were renor- 
malized to the data in order to remove the energy-de- 
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Fig. 12. Main picture: angular distributions for incident pho- 
ton energies from 152 MeV (bottom curve) to 264 MeV (top 
curve, same energy bins as in Fig. I10|l compared to PWIA 
results using the fitted form factor. Absolute scales of model 
results normalized to data (see text). The insert shows on a 
larger scale the low-energy results (144 MeV - 168 MeV). 



pendent FSI effects. This PWIA must describe by con- 
struction the angular distributions on average. However, 
it actually agrees almost perfectly with the shapes of all 
individual distributions, with very different relations be- 
tween pion angles and nuclear momentum transfers. This 
means that an FSI effect would be needed, which over a 
range of incident photon energy of more than 100 MeV 
has exactly the same angular and momentum-transfer de- 
pendence as a change of the form factor from an rms ra- 
dius of 2.3 fm to 2.66 fm. Although this does not seem 
to be a likely scenario, reasonably sophisticated modelling 
of the FSI effects is needed before a final conclusion can 
be drawn. However, results for a similar analysis of coher- 
ent photoproduction off carbon, calcium, and lead nuclei 
point to a small influence of FSI on the extracted radii for 
light nuclei. Fully taking into account the FSI effects [5] 
lowered the extracted value of the mass radius for 208 Pb 
by 5.8%, for 40 Ca by 2.2%, but for the lighter 12 C only by 
0.9%, while the observed difference between charge and 
mass radius for 7 Li is on the 10% level. 
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5.2 Coherent ^-photoproduction 

The total cross sections extracted for the two 77-decay 
channels, summarized in Fig. 1131 are nicely consistent. 
They show a much smoother rise at production threshold 
than the 3 He data (cf. Fig.QJ. For a quantitative analysis 
their average is compared to PWIA modelling, based on 
Eq. [JJ in Fig. [TU As discussed in Sec. [5] the situation is 
much different to pion production since for ^-production 
the cross section is dominated by the contribution of the 
odd lj>3/2 proton, which is only a small correction in the 
7T° case. 




660 



ETMeV] 



Fig. 13. Comparison of the total cross section for coherent re- 
production from the two-photon and six-photon decay of the r\. 
The vertical dotted lines indicate coherent and breakup thresh- 
olds. 



Results from PWIA, using the charge form factor or 
the mass form factor fitted to the pion data (dashed, 
respectively dotted curves in Fig. are similar. This 
is simply so because, for 77 production, large momentum 
transfers dominate where the two form factors agree. The 
relative contribution of the incoherent excitation of the 
478-keV state is significant in PWIA. Also this had to be 
expected because for ^-production there is no piece with 
an A 2 -term like Eq. 0] for coherent pion production and 
the elastic and inelastic form factors are similar for large 
momentum transfers. The systematic uncertainty of the 
PWIA results is larger than in the pion case because the 
cross section is dominated by these less well-established 
contributions. However, altogether the comparison of the 
energy dependence of the measured total cross section and 
the PWIA results in Fig. [HI shows no threshold enhance- 
ment above phase-space behavior, and thus no indication 
for the formation of a quasi-bound state. The situation 
is thus much different from the 3 He case discussed in the 
introduction which, apart from the incoherent excitation, 
has similar systematic uncertainties in PWIA. Compar- 
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Fig. 14. Comparison of the average of the experimental two- 
photon and six-photon cross sections to the PWIA results. 
Dashed (dotted) curves: coherent contribution (Eq. [7] without 
F^^-term) based on charge form factor (mass form factor fit- 
ted to pion production). Solid (dash-dotted) curves: sum of co- 
herent and incoherent contribution (see Eq.[7J for charge (mass 
form factor). The insert shows the ratio of measured cross sec- 
tion and PWIA results, open symbols only coherent part, filled 
symbols sum of coherent and incoherent contributions). 
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Fig. 15. Comparison of measured angular distributions for 
7+ 7 Li— > 7 Li+?7 (red) dots to PWIA results, solid lines: full 
PWIA, dashed lines: only coherent part. 



ison of the two results highlights the special role of the 
rj — 3 He system. 

Also the results for the angular distributions, summa- 
rized in Fig. 1151 are consistent with this interpretation. 
They agree better with the momentum-transfer depen- 
dence of the form factor than in the 3 He case and show no 
tendency towards isotropic behavior close to threshold. 
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6 Summary and Conclusions 

Precise data have been measured for coherent photopro- 
duction of 7r°-mesons off 7 Li nuclei and coherent photopro- 
duction of 77-mesons off the same nucleus has been identi- 
fied for the first time. The experimental results for the pion 
production are quite well reproduced at low incident pho- 
ton energies by a PWIA dominated by the spin/isospin- 
independent part of the elementary production amplitude. 
The spin- flip amplitude from the unpaired 1^3/2 proton is 
considered for the leading M\ + multipole and the correc- 
tions applied for the incoherent excitation of the 478-keV 
nuclear state in 7 Li are insignificant. This model repro- 
duces quite well total cross sections and angular distribu- 
tions at incident photon energies below 225 MeV, indicat- 
ing that distortion effects from final-state interactions of 
the pion are small in this energy range. 

After an adjustment of the nuclear form factor, which 
corresponds to a change in the harmonic double-well pa- 
rameterization from an rms radius of 2.3 fm (reported for 
the charge form factor derived from electron scattering 
data) to 2.66 fm, the shape of angular distributions in 
this energy range is excellently reproduced. Exploiting the 
possible uncertainties due to approximations, in particu- 
lar in the spin-flip term of the PWIA, we find reasonable 
agreement between data and PWIA for rms radii down 
to 2.5 fm, which are still larger than previously reported 
charge radii and also larger than predictions for the mass 
radius, which are around 2.35 fm (see e.g. [7,33 ). DWIA 
calculations with careful treatment of possible FSI effects 
are needed for further analysis of this discrepancy. 

Coherent photoproduction of ?7-mesons is quite diffi- 
cult to measure and so far only results for the deuteron 
[231124] and 3 He [91HU] had been reported. This experiment 
extended the mass range to 7 Li by measuring total cross 
sections on the level below 20 nb. The results, also for the 
angular distributions, are in good agreement with PWIA 
expectations and do not show an unexplained threshold 
enhancement as in the 3 He case, underlining the special 
role of 3 He as a candidate for ^-mesic states. 
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